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Type la supemovae are key tools for measuring distances on a cosmic scale. 
They are generally thought to be the thermonuclear explosion of an accreting 
white dwarf in a close binary system. The nature of the mass donor is still 
uncertain. In the single-degenerate model it is a main-sequence star or an 
evolved star, whereas in the double-degenerate model it is another white dwarf. 
We show that the velocity structure of absorbing material along the line of 
sight to 35 type la supemovae tends to be blueshifted. These structures are 
likely signatures of gas outflows from the supernova progenitor systems. Thus 
many type la supemovae in nearby spiral galaxies may originate in single- 
degenerate systems. 

Type la supemovae (SNe la) have large and calibratable luminosities, making them good 
tools for measuring distances on a cosmic scale to gauge the geometry and evolution of the 
Universe (i, 2). Understanding the nature of the progenitor system is important, as progenitor 
evolution or a changing mix of different progenitors may bias cosmological inferences. The 
consensus view of SNe la is that mass transfer onto a massive carbon-oxygen white-dwarf (WD) 
star in a close binary leads to a thermonuclear explosion, as the mass of the WD approaches 
the critical Chandrasekhar mass limit (5). In the single-degenerate (SD) model the mass donor 
is either a main-sequence star or an evolved subgiant or giant star, whereas in the competing 
double-degenerate (DD) model it is another WD (4). 

In the SD scenario, nonaccreted material blown away from the system before the explosion 
should remain as circumstellar matter (CSM) (5); thus, detection of CSM in SN la spectra would 
lend support to the SD model. Patat et al. reported such a detection based on time-variable 
absorption features of the Na I doublet (Na I Di and D2; A^i = 5896 A and Xd^ = 5890 A 
rest-frame wavelengths) in optical spectra of SN 2006X (6). These authors suggest that CSM 
may be common to all SNe la, though variation in its detectability can exist due to viewing 
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angle effects. Although multi-epoch high-resolution spectra of SN 2000cx (7) and SN 2007af 
(8) exhibited no absorption features or no time variation in them, additional detections were 
recently reported for three other events: SNe 20071e, 1999cl, and 2006dd (9-11). While it has 
been suggested that Na D absorption could not have been caused by CSM due to ionization 
considerations (72), detailed arguments to the contrary have been presented (9, 13). 

Absorption features from unrelated intervening gas clouds are expected to have random 
velocity offsets with respect to the SN, whereas absorption due to winds blown by the progenitor 
system is expected to be always blueshifted with respect to the SN, as the source of light is 
behind the outflowing absorbing material. We thus searched for this signature in a large sample 
of high-resolution, single-epoch observations of SNe la. Our sample consisted of spectra of 35 
SNe la and 1 1 core-collapse (CC) SNe obtained using the Keck HIRES and Magellan MIKE 
spectrometers. We also studied previously published spectra of six SNe la (6, 14-17) and seven 
CC SNe (16-22). A bias may exist if other spectra were obtained but not published. Therefore, 
we report results separately with and without these historical SNe (hereafter "the extended 
sample", for which results are given in square brackets). Detailed analyses of a few events from 
our sample have been pubUshed: SN 2006X (6), SN 2007af (7), SN 20071e (9), and SN 2008D 
(23). 

The HIRES observations used setups with spectral resolution (R = X/6X) of ~ 50,000- 
52,000, giving full width at half-maximum (FWHM) intensity of ~ 0.115-0.12 A (or ~ 5.5-6 
km s~^) in the vicinity of the Na I D lines. The MIKE observations were obtained with R ^ 
30,000-36,000, giving FWHM « 0.165-0.2 A (or - 8.4-10 km s"^) near Na I D. We reduced 
the HIRES spectra using the MAKEE (MAuna Kea Echelle Extraction) pipeline (24). The 
MIKE data were reduced using the latest version of the MIKE pipeline (25). 

We obtained the redshift and morphological classification of the host galaxies of our SN 
sample from the NASA/IR\C Extragalactic Database (NED; 26), the SIMBAD astronomical 
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database (27), or using images from the Digital Sky Survey (DSS; 28). Our sample host-galaxy 
redshifts z are between 0.0019 [0.0015] and 0.06. Tables SI and S2 list our complete SN sample 
with host-galaxy properties. 

The redshifts of the host galaxies are sufficiently large to allow us to differentiate between 
absorption resulting from material in the Milky Way (MW; i.e., at 2; ~ 0; used as a control 
sample) and from material in the host galaxies of the SNe. Out of the 35 [41] events in our 
SN la sample, 22 [28] events exhibit absorption features consistent with the object's host-galaxy 
redshift and 13 events do not [supporting online material (SOM) text SI]. In the CC SN sample, 
9 [16] events exhibit absorption compatible with the SN host-galaxy redshifts and 2 events do 
not (SOM text SI). Of a total of 46 [59] events (SNe la and CC SNe), 42 [51] exhibit Na I D 
absorption due to material in the Milky Way whereas 4 [8] do not. 

We normalized the spectra across the continuum and translated wavelengths into velocities 
relative to the wavelength of the minimum of the most prominent absorption feature of both of 
the D lines (SOM text S2), both for the features in the host galaxies and in the Milky Way when 
detected, using the Doppler shift, 

^;/c^(A-Ao)/Ao, (1) 

where c is the speed of light in vacuum, Aq is the wavelength of the zero-velocity component, 
and V is the velocity of the component that was Doppler shifted from Aq to A. We categorized 
the absorption features exhibited in our spectra into three classes of structures as follows (see 
Fig. 1 for graphic examples): 

(i) Blueshifted: One prominent absorption feature with weaker features at shorter wave- 
lengths with respect to it. 

(ii) Redshifted: One prominent absorption feature with weaker features at longer wavelengths 
with respect to it. 
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Fig. 1: Examples of spectra classified as blueshifted (top left), redshifted (top right), and single 
(bottom left) or symmetric (bottom right). Flux is normalized and velocities are given in km s~^ 
relative to the prominent absorption-feature minimum (relative to D2 in blue and to Di in red). 

(iii) Single/Symmetric: A single absorption feature, or several features with both blue and 
redshifted structures of similar magnitude. 



Classification results are presented in Table 1, Figure 2, and Figure SI. 



We used the Galactic Na I D absorption along the line of sight to our SNe as a control 
sample, extending it with Galactic absorption along the line of sight toward quasi-stellar objects 
(QSOs) published by Ben Bekhti et al. (29). To imitate being observed from outside of the 
Milky Way we inverted the velocities of the Galactic features. Figures S2-S5 show the spectra 
of the SN la sample and Figure S6 shows the spectra of the CC SN sample. 

The absorption features we observe are narrow, with typical separations and velocity dis- 
persions of a few tens of km s^^, and are distinct from features due to material ejected in 
the explosion itself that are wide (~ 10^ km s^^), reflecting the expansion velocity of the SN 
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Fig. 2: Pie-chart visualization of the results given in Table 1. The results for the extended 
samples are shown in Fig. SI. 

ejecta. It is apparent from Table 1, Figure 2, and Figure SI that the SN la sample displays a 
strong preference for blueshifted structures, whereas the Milky Way sample shows no signifi- 
cant preference. 

To test whether our SN la sample could be a random draw out of a uniform distribution, 
we calculated the probability of observing a set of 22 [28] spectra exhibiting the preference we 
detected in our sample, or a set even more extreme (i.e., 12 [16] or more blueshifted spectra 
and 5 [6] or fewer redshifted spectra); the result was low 2.23% [0.54%]. Moreover, a K-S test 
rejects the null hypothesis that the SN la and MW samples are from the same distribution at a 
1% significance level. 

Out of the 17 SNe la that were classified as blueshifted or redshifted, 2 events (1 1.8%; SNe 
2006X and 2008fp) exhibit saturated features. Out of the 7 CC SNe classified as blueshifted 
or redshifted, 4 events (57.1%; SNe 2006bp, 2008cg,x 2008D, and 2008 J) exhibit saturated 
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Table 1: Classification of absorption features. 



Sample 


Blueshifted 


Redshifted 


Single/Symmetric 


Total 


SNe la 


12 [16] 


5 [6] 


5 [6] 


22 [28] 


CC SNe 


4 [8] 


3 [5] 


2 [3] 


9 [16] 


MW (SNe) 


12 [13] 


13 [16] 


17 [22] 


42 [51] 


MW (QSO) 


10 


13 


6 


29 


MW (QSO+SNe) 


22 [23] 


26 [29] 


23 [28] 


71 [80] 



features. Three SNe la in our sample occurred in elliptical galaxies (SNe 2006ct, 2006eu, and 
2007on; see Table SI), and their spectra do not exhibit Na I D absorption features. Eight of our 
SNe la occurred in lenticular galaxies, of which two exhibit absorption features (one blueshifted 
and one single) and 6 do not. Therefore, our sample provides a good representation of SNe la 
in nearby (z < 0.06) spiral galaxies. 

We calculated the velocities and column density (N) of the absorbing material by fitting 
them with Voigt profiles useing VPFIT [(30); see Figure 3]. SN la features have a log(A'^) similar 
to the that of the MW features, lower on average than those in CC SNe. The relative velocities of 
the SN la features range is between —150.8 to +139.4 km s~^ (with an additional single feature 
at +193.7 km s~^), whereas the CC SN relative velocities range is between —87.2 to +116.6 
km s~^, and the MW velocity range is between —79.5 to +116.9 km s~^. The cumulative 
distribution of the relative number of features as a function of the velocity for each sample is 
given in Fig. 4. 

The narrow absorption features we measured could possibly be produced by, 

(i) intervening clouds of gas in the host galaxy; 

(ii) a systemic wind blown by the host galaxy; 
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Fig. 3: A plot of log(A^) of the absorbing material calculated using VPFIT as a function of 
velocity relative to the zero-velocity component, for the blueshifted and redshifted components 
observed in our sample. Zero-velocity components are omitted. The solid lines are the log(A^) 
averages for the different samples: SNe la (blue), log(A^)avg = 11.67 ± 0.054; CC SNe (red), 



log(iV) 



avg 



12.07 ± 0.075; and MW (green), log(A^) 



are given in cm 
lines. 



11.59 ± 0.058. Column densities 
Dashed lines mark the uncertainty in the mean of the corresponding solid 
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Fig. 4: A plot of the cumulative density function of our three samples (SNe la, CC SNe, and 
MW). The zero- velocity components are omitted in order to minimize the effect of the different 
sample sizes. The color scheme is the same as that used in Figure 3. The overabundance of 
SN la high-velocity bluesifted absorption components (—150 km s^^ < Ujci < —100 km s^^) 
in comparison with the CC SN sample is questionable, as a K-S test cannot reject (at a 5% 
significance level) the null hypothesis that both samples are from the same parent distribution. 
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(iii) nonrandomly moving interstellar matter (ISM) - a relic from a star-forming event; and 

(iv) material ejected from the progenitor system before its explosion. 

In the first case, intervening clouds of gas in galaxies may have a small random velocity relative 
to the local rotation velocity; therefore, the resulting absorption features are equally likely to 
be either blueshifted or redshifted relative to the SN, in contrast to the preference exhibited 
in our SN la sample (shown to be inconsistent with a random draw from a uniform parent 
distribution). The structure of the absorption features observed in the Milky Way is due to the 
movement of ISM gas clouds, and our results are in good agreement with a uniform distribution 
from randomly moving clouds (see Table 1, Fig. 2, and Fig. SI). The CC SN sample is 
unfortunately quite small, and is consistent with either a uniform distribution or one similar to 
that of the SN la sample. 

If the second explanation was viable, we would expect to see a signature of this wind in the 
Milky Way sample, as a prevalence of blueshifted absorption (i.e., outflowing). However, the 
MW could be a special case. We would also expect to see a similar, or even greater, effect in 
our CC SN sample, but this is not seen. Nevertheless, the CC SN sample is too small to reach 
statistical significance. Furthermore, if this model was correct, blueshift due to galactic outflows 
would be common, and imply a large velocity span, at small host inclination (i.e., hosts closer 
to being face on). Figure S7 shows that the absorption features of SN la events that occurred in 
spiral hosts (see Figure S9-S11 for host images) show no preference to low-inclination hosts, 
nor is there a correlation between the velocity span and the inclination. 

In addition, galactic winds are observed mainly in galaxies with a star-formation-rate (SFR) 
surface density threshold above ~ O.IMq yr~^ kpc~^, and are rare in local ordinary disk galax- 
ies [(31) and references therein]. Moreover, galactic winds are thought to be expelled into the 
halo and then cool and fall back onto the disk [the "galactic fountain" process; (32-33)]. Be- 
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cause of its low first ionization potential (5.1 eV) sodium is a good tracer for neutral gas; thus, 
Na 1 D absorption is expected to be observed in the later stage of the galactic fountain process 
(i.e., the inflow) when the gas is mostly neutral, and less during the earlier outflow stage when 
the gas is expected to be mostly ionized (34-35), resulting in a redshifted preference. Absorp- 
tion from gas in galactic halos has larger equivalent widths and column densities of Ca 11 H & 
K, A^(Ca H) relative to Na I D, while ISM and CSM have A^(Na l)/M(Ca II) of order unity [(36 
and references therein]. We have spectral coverage of Ca II H & K for 12 out of the 22 SN la 
events that exhibit Na 1 D absorption (seven of which are blueshifted). All these events exhibit 
A/^(Na I)/AA(Ca II) < 1. This implies that the absorbing material is not in the galactic halo. 

It is plausible that most of the SNe in our SN la sample are "prompt" events, arising from 
progenitor systems with short (< 500 Myr) delay times from star formation to explosion (37- 
38). If this is the case, the massive stars that were formed at the same time as the progenitor 
system could have driven a wind and blown a cavity in the ISM, such that the remaining gas 
clouds along the line of sight would be blueshifted from our perspective. However, even in spiral 
galaxies, SNe la do not follow massive stars as traced by H II regions (39), and so are unlikely to 
reside close to recent star-formation regions. This process might be evident in our observations 
of CC SNe which show an excess of strong (and even saturated) absorption compared to the 
SNe la. 

Hence, we conclude that the blueshifted preference is due to an inherent feature of the SNe 
themselves. The low velocities and narrowness of the features suggest that these absorption 
features are due to CSM that was ejected by the progenitor system prior to the SN explosion. 
Features with relatively low velocities (i.e., of order a few 10 km s~^) may be attributed to a 
wind blown by a red giant companion, but those with higher velocities are incompatible with 
such winds. 

Patat et al. (6) proposed that the observed features in spectra of SN 2006X arise from shells 



11 



of successive recurrent nova eruptions that swept up the stellar wind of the companion [SOM 
text S6 of (6)]. This scenario has gained support from a study of the Na 1 D, Ca 11 H&K, and 
K I blueshifted absorption lines observed in spectra of RS Ophiuchi, a symbiotic recurrent nova 
(40). These features are similar to the features exhibited in our SN la sample (SOM text S3), 
i.e., at least one type of SD system is shown to be able to produce such features. 

Our findings are consistent with the SD model (41-42) as the progenitor system of a certain 
fraction of SNe la in nearby spiral galaxies (43). But how large is this fraction? The detection of 
the absorption features might be angle dependent; thus, we can only estimate a lower limit. If we 
assume that we should observe the same amount of blueshifted and redshifted components due 
to ISM, we can attribute any blueshifted excess to CSM, implying at lower limit of 25% [26%] 
of SNe la events that occur in spiral galaxies. Alternatively, assuming the MW blueshifted 
fraction represents galactic ISM, we get an estimate of at least 20% [24%] of SNe la in spiral 
hosts. 
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Supporting Online Material 

SI. 3cr Upper Limits for Detection of Absorption Features when Features are Not Appar- 
ent 

Thirteen of the SNe la and two CC SNe in our sample do not exhibit Na I D absorption 
features. Using Equation 4 of Leonard & Filippenko (45), we calculated Wx{3a), the upper 
detection limit of Na I D absorption for these spectra: 



where AA is the spectral resolution, A/ is the root-mean square (rms) fluctuation of the flux 
about the normalized continuum level, Wune is the spread of the D2 or Di lines (in our sample, 
we see that the Na I D features are usually within a spread of ~ 200 km s"^ or ~ 4 A), and 
B is the number of bins per resolution element. We used AA = 0.115 A, A/ = (S/N)~^, 
Inline = 4 A, and S = 5; we calculated the signal-to-noise ratio (S/N) over a range of 2 A of 
continuum between the D2 and Di lines using the IRAF splot tool. Table S3 lists the values of 
S/N and Wx{3a) for the nondetection events. 

Among the spectra that exhibit Na I D absorption features, those of SNe 2006E, 2007 af, and 
2007rw have the weakest ones (see Fig. S2, S3, and S6). Fig. S12 shows the Gaussian fits to 
the D2 lines of these events using the parameters given in Table S4. It is apparent from Table S3 
that in the events whose S/N is above 10, features like those exhibited by SNe 2006E or 2007af 
should have been apparent; hence, we conclude that any Na I D absorption must be weaker than 
any of the detected events. The line in the SN 2007rw spectrum has an equivalent width that is 
less than Wx{3a) (=54.9 mA). 




(2) 
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52. The Zero- Velocity Component 

Ideally, we would like to know the progenitor system's velocity in order to use it to define 
our zero- velocity frame. However, this velocity is not easily determined. In some of the cases 
(SNe 2006cm, 2006X, 2007af, 20071e, 2007sr, 2008C, 2009ds, and SNF 20080612-003), a 
distinct Ha emission line (from superposed H 11 regions) can be seen in the spectra, enabling 
us to determine the local frame of reference (LFR) of the galaxy along the line of sight. These 
features cannot be attributed to the SNe la explosions themselves as they are narrow emission 
lines. In most of the cases the progenitor system should have a velocity close to that of the LFR. 
When present, the Ha emission, and thus most probably the LFR, coincides (within a few km 
s"^ to a few 10 km s~^, in all but one case) with the velocity of the most prominent Na I D 
absorption feature, and blueshifted components extend out to 150 km s~^. 

In cases when Ha is not detected, we calculate the expected location of the Na I D lines 
using the cataloged galaxy redshifts from NED (26), and find these to coincide with the most 
prominent Na I D absorption to within a few tens of km s~^; weaker features extend up to 200 
km s^ (46). Though we cannot determine the absolute velocity of the blueshifted absorption 
with respect to the progenitor system in each individual case, we can say with high confidence 
that these velocities span a few tens of km to 200 km s~^. Note that the strong preference 
we see in SNe la for dominant strong features with bluer weaker features does not depend on 
the location of the velocity zero point. This preference is naturally explained by progenitor 
outflows and inconsistent with random ISM clouds. 

53. Velocity of Recurrent Novae Shells 

As a first-order approximation, Patat et al. (40) show that a shell formed by nova ejecta 
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sweeping up the stellar wind of the companion star expands at a velocity of roughly 

u = ——^ , (3) 

1 + A* 

where vq is the scale velocity of the nova ejecta (p oc exp{—v/vo)), Uy, is the velocity of the 
companion's wind, and fj, = Me/M^ is the mass ratio between the nova ejecta mass and the 
swept-up mass. They show that for typical values of these parameters the velocity of such a 
shell will be ~ 50 km relative to the systemic velocity of the recurrent nova system. 

When the zero velocity is taken to be the minimum of the deepest absorption feature, as 
is done in our analysis, the bluest feature of RS Oph has a velocity of ~ 70 km s~^ and the 
features show resemblance to features in our SN la sample [see fig. 10 of (40)]. 

Nevertheless, some of our SN la events exhibit features with higher velocities. The con- 
sistency of these features with this first-order approximation needs to be checked. Using the 
same parameters used for the slower feature (Me = 10~^ Mq, vq — 9.1 x 10^ cm s~^, — 20 
km s"^), and assuming m = 150 km s"^ we get n = 0.0778, or M«, 1.3 x 10~^ Mq. This 
value implies a red giant mass loss of Mrg ~ 1.3 x 10"^/ At w 6.5 x 10~^ Mq yr~^ (the latter 
assuming At — 20 yr) which is a reasonable value (47). 
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Table S 1 : SN la sample host-galaxy morphology classification 



Absorption classification 


Object name 


Host galaxy 


Type 




SN 2006X 


NGC 4321 


Sbc 




SN 2007fb 


UGC 12859 


Sbc 




SN 2007kk 


UGC 2828 


Sbc 




SN 20071e 


NGC 7721 


Sc 




SN 2008C 


UGC 3611 


so/a 


Blueshifted 


SN 2008dt 


NGC 6261 


so/a 


SN 2008ec 


NGC 7469 


Sa 




SN 2008fp 


ESQ 428-G14 


SO 




SNF200806 12-003 


2MASX J16152860-1-1913344 


Sc* 




SN 2009ds 


NGC 3905 


Sc 




SN 2009ig 


NGC 1015 


Sa 




SN 2009iw 


IC 2160 


Sbc 




SN 2007af 


NGC 5584 


Scd 




SN 2007sa 


NGC 3499 


so/a 


Redshifted 


SN 20091e 


ESQ 478-G006 


Sbc 




SN 2009mz 


NGC 5426 


Sc pec 




SN 2010A 


UGC 2019 


Sab 




SN 2006cm 


UGC 11723 


Sb edge-on 




SN 2006E 


NGC 5338 


SO 


Single/ symmetric 


SN 2007fs 


ESQ 601-G5 


Sb 




SN 2007sr 


NGC 4038 


Smpec 




SN 2010ev 


NGC 3244 


Scd 




SN 2006ct 


2MASX J 12095669-1-4705461 


E* 




SN 2006eu 


MCG -1-08-36-16 


E 




SN 2007gi 


ESQ 298-28 


Sbc 




SN 2007hj 


NGC 7461 


SO 




SN 2007on 


NGC 1404 


El 


None 


SNF200805 14-002 


UGC 8472 


SO 


SN 2008dh 


2MASX J0035 1050-1-23 15 184 


Sb* 




SN 2008ee 


NGC 307 


SO 




SN 2008ge 


NGC 1527 


SO 




SN 2008hv 


NGC 2765 


so 




SN 2008ia 


ESQ 125-G006 


so 




SN 2009ev 


NGC 5026 


so/a 




SN 2009nr 


UGC 8255 


Scd 



*Morphology classified using images from DSS (Digital Sky Survey) (28). 



21 



Table S2: CC SN sample host-galaxy morphology classification 



Absorption classification 


Object name 


Host galaxy 


Type 


Diucsniiicii 


SN 2006bp 


NGC 3953 


Sbc 


SN 2006CU 


UGC 9530 


S 


SN 2008ax 


NGC 4490 


Sd pec 


SN 2008cg 


FGC 1965 


Scd 


Redshifted 


SN 2006be 


IC 4582 


S 


SN 2008D 


NGC 2770 


Sc 


SN 2008J 


MCG -02-7-33 


Sbc 


Single/Symmetric 


SN 2003gd 


NGC 628 


Sc 


SN 2007rw 


UGC 7798 


Im 


None 


SN 2006ca 


UGC 11214 


Scd 


SN 2007pk 


NGC 579 


Scd 



Table S3: S/N and Wx{3a) for spectra with no apparent Na I D absorption features 



Object name 


Type 


S/N* 


Wx{3(t) (mAy* 


SN 2006ct 




6.47 


144 


SN 2006eu 




11.11 


84 


SN 2007gj 




12.52 


74 


SN 2007hj 




7.22 


129 


SN 2007on 




39.95 


23 


SNF200805 14-002 




26.6 


35 


SN 2008dh 


la 


7.52 


124 


SN 2008ee 




23.41 


40 


SN 2008ge 




111.27 


15.5 


SN 2008hv 




85.94 


20 


SN 2008ia 




55.03 


31.4 


SN 2009ev 




76 


16 


SN 2009nr 




108.53 


15.9 


SN 2006ca 


II 


22.30 


42 


SN 2007pk 


Iln-pec 


6.47 


144 



*Calculated over a range of 2 A of continuum between the D2 and Di lines using the IRAF splot 
tool. 

**For each one of the Na I D lines. 
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Table S4: Gaussian fit parameters for SN 2006E and SN 2007 afPa features 



Object name 


Ao(A) 


FWHM (mA) 


W(mA) 


S/N 


SN 2006E 


5906.7 


192 


92 


12.21* 


SN 2007af 


5921.7 
5922.1 


286 
168 


127.5 
39.3 


47.72** 


SN 2007rw 


5941.54 


115 


21 


16.57*** 



*Calculated over the range 5908-5910 A. 
**Calculated over the range 5924-5926 A. 
***Calculated over the range 5943-5945 A. 
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Figure S2: SN la spectra that exhibit blueshifted structures. The blue/red color scheme is as 
described in the caption of Fig. 1. Two additional spectra of SN 1991T and SN 1998es are 
given in Fig. S5 of (6). Spectra of SN 2006X and SN 20071e from additional epochs can be 
found in Patat et al. (2007) {6; structure of blueshifted components seen better) and Simon et 
al. (2009) (9), respectively. 24 
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Figure S3: SN la redshifted structure spectra. Color scheme identical to that described in the 
caption of Fig. 1. Spectra of SN 2007 af from additional epochs can be found in Simon et al. 
(2007) (8). 
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Figure S4: SN la spectra that exhibit single/symmetric absorption structures. The blue/red color 
scheme is as in Fig. 1. 
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Figure S5: SN la spectra that exhibit no Na I D absorption. The color scheme is as described in 
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Figure S6: Spectra of CC SN sample. Upper four panels are of events exhibiting blueshifted 
structures. Middle three panels are of events exhibiting redshifted structures. The bottom four 
panels are of events classified as single/symmetric features (SN 2003gd and SN 2007rw) and 
those that do not exhibit Na I D absorption features (SN 2006ca and SN 2007pk). The color 
scheme is as described in Fig. 1. Analysis of SN 2008D has been reported by Soderberg et al. 
(2008) (2i). 
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Figure S7: The velocity span of the SN la absorption features as a function of host galaxy 
inclination for events with spiral hosts (inclination 0° is face on and 90° is edge on). Inclinations 
are calculated using the major and minor axis given in NED (26), cos{i) = h/a (where i is the 
inclination, a the major axis, and b is the minor axis). Blue lines represent the events classified 
as blueshifted. Red lines represent redshifted events and green lines represent symmetric 
events. Green squares represent events with a single feature and black circles represent events 
with no sodium absorption detection. Galactic outflows should exhibit blueshifted absorption 
with a larger velocity span at lower host galaxy inclination. This is not observed in our SN la 
sample. There seems to be no correlation between the host galaxy inclination and the velocities 
of the absorption features as expected from CSM. 
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Figure S8: Images of the spiral host galaxies of blueshifted SNe la. The SN names and the host 
inclinations are given above each frame. Dark background images were taken from the ESO 
DSS archive (http://archive.eso.org/dss/dss). Light background images were taken using the 
SDSS DR8 finding chart tool (http://skyserv^ftdss3.org/dr8/en/tools/chart/chart.asp). Scales 
are indicated at the upper left comer of each image. Cross-hairs mark the location of the SNe. 






Figure S9: Same as S8 for the redshifted SNe la. 
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SN 2007gj, i = 62° SN 2008dh, i = 71 




SN 2009ev, i = 51° SN 2009nr, i = 40' 




Figure S 1 1 : Same as S8 for the SNe la with no detection of sodium absorption. The inclination 
of the host galaxy of SN 2008dh was calculated by fitting an ellipse to the host and calculating 
the ratio between its minor to it's major axis. 
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Figure S12: Gaussian fits of the Na I D2 absorption features exhibited in spectra of SN 2006E 
(upper), SN 2007af (middle), and SN 2007rw (lower). Parameters of the fits are give in Table 
S4. Blue dots are the observed data. In the upper figure the black line is the Gaussian fit to the 
single feature. In the lower figure the red/green line is the fit to the left/right absorption feature, 
and the black line is the fit to both features combined. The vertical lines mark the center of the 
Gaussian fits. 34 



